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ABSTRACT 

^ \ We study the large-scale stellar and gaseous kinematics of the E4 galaxy NGC 2974, based 

on panoramic integral-field data obtained with SAURON . We quantify the velocity fields with 
Fourier methods (kinemetry), and show that the large-scale kinematics is largely consistent 
with axisymmetry. We construct general axisymmetric dynamical models for the stellar mo- 
tions using Schwarzschild's orbit-superposition method, and compare the inferred inclination 
\ and mass-to-light ratio with the values obtained by modelling the gas kinematics. Both ap- 

. proaches give consistent results. However we find that the stellar models provide fairly weak 

constraints on the inclination. The intrinsic orbital distribution of NGC 2974, which we in- 
fer from our model, is characterised by a large-scale stellar component of high angular mo- 
mentum. We create semi-analytic test models, resembling NGC 2974, to study the ability of 
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O ' Schwarzschild's modelling technique to recover the given input parameters (mass-to-light ra- 

|H , tio and inclination) and the distribution function. We also test the influence of a limited spatial 

c/3 ■ coverage on the recovery of the distribution function (i.e. the orbital structure). We find that 

^ ' the models can accurately recover the input mass-to-light ratio, but we confirm that even with 

perfect input kinematics the inclination is only marginally constrained. This suggests a possi- 
ble degeneracy in the determination of the inclination, but further investigations are needed to 
rS \ clarify this issue. For a given potential, we find that the analytic distribution function of our 

^ . test model is well recovered by the three-integral model within the spatial region constrained 

by integral-field kinematics. 

Key words: galaxies: elliptical and lenticular - galaxies: kinematics and dynamics - galaxies: 
structure, galaxies: individual, NGC 2974 



1 INTRODUCTION 

The internal dynamical structure of galaxies retains evidence of 
their evolution. The internal dynamics, however, can only be inter- 
preted through a combination of observational and theoretical ef- 
forts. From a theoretical point of view, one wants to know how the 
stars are distributed in space and what velocities they have. From 
the observational point of view, one wants to determine the intrin- 
sic structure of the observed galaxies. The goals of both approaches 
are equivalent, and consist of the recovery of the phase-space den- 
sity, or distribution function (DF) of galaxies, which uniquely spec- 
ifies their properties. An insight into the DF is possible by the 
construction of dynamical models which are constrained by ob- 
servations. There are several modelling methods established in 
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the literature, of which Schwa rzschild' s orbit-superposition method 
is perhaps the most elegant iSchwarzschil rf' 197^. Tisl. In the 
past few ye ars it has been applied successfully t o a number 
of galaxies I van der Marel et al. 1998; Cretton & v an den BoschI 
.1999; Cappellari et aL.2002: .Gebhardt et aL 2003) : but recent ob- 
servational advances in spectroscopy with integral-field units of- 
fer for the first time full two-dimensional constraints on these dy- 
namical models ( Verolme et al. 2002; Copin, Cretton & Emsellem 
2004). 

This paper presents a case study of the early-type galaxy 
NGC 2974. It is one of the few elliptical galaxies known to 
contain an extended disc of neutral hydrogen in regular ro- 
tation iKimetal. 198"3). It also hosts extended Ha emission 
(Busonetal. 1993; Plana etal. 1998), and belongs to the 'rapid 
rotators' liBendeik.l988i) . The total absolute magnitude of Mb = 
—20.32 puts NGC 2974 near the transition between giant ellip- 
ticals and the lower-luminosity objects which often show pho- 
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tometric a nd kinematic evide n ce for a significant dis c compo- 
nent (e.g. lRix&Whitj|l992h . lEmsellem et"ai] J2003l hereafter 
EGF03), combining WFPC2 imaging with TIGER integral-field 
spectroscopy of the central few arcseconds, discovered a spiral 
structure in the Ha emission in the inner few arcseconds, and con- 
cluded that the galaxy contains a inner stellar bar. The general prop- 
erties of NGC 2974 are listed in TableQ 

The availability of both stellar and gaseous kinematics makes 
NGC 2974 a very interesting ca s e for detailed dynamical mod- 
elling. |c inzano & van der Marell ^993) made dynamical Jeans 
models of the gaseous and stellar components additionally intro- 
ducing a stellar disc in order to fit their long-slit data along three 
position angles. They found that the stellar and gaseous discs were 
kinematically aligned and the inclination of both discs was consis- 
tent with 60° . This prompted them to suggest a common evolution, 
where the gas could be ionised by the stars in the stellar disc. Us- 
ing more sophisticated two-integral axisymmetric models, which 
assume the DP depends only on the two classical integrals of mo- 
tion, the energy E and the angular momentum with respect to the 
symmetry axis L^, EGF03 were able to reproduce all features of 
ICinzano & van d er Marel' ( 1994) data as well as their integral-field 
TIGER data (covering the inner 4"). The models of EGF03 did not 
require a thin stellar disc to fit the data. 

In this study we construct axisymmetric models for NGC 2974 
based on Schwarzschild's orbit superposition method. This method 
allows the DF to depend on all three isolating integrals of motion. 
All previous studies with three-integral models concentrated on the 
determination of the mass-to-light ratio T, and mass of the central 
black hole, A/bh. Based on the observed stellar velocity disper- 
sion, the A/bh — cr relation (e.g. Tremaine et al. 2002) predicts a 
central black hole mass o f 2.5 x Mp,, which at the distance 
of NGC 2974 (21.48 Mpc. lTonrv et all200ll) has a radius of influ- 
ence of 0'.'2. Our obser vations of NGC 297 4, with the integral-field 
spectrograph SAURON jBaconetall200lh . do not have the neces- 
sary resolution to probe the sphere of influence of the central black 
hole. The dynamical models presented here are therefore aimed at 
determination of the T, the inclination, i, and the internal orbital 
structure. The stellar and gaseous kinematics also provide indepen- 
dent estimates of T and i, which can be used to cross-validate the 
results from the two approaches. 

The results of the dynamical modelling are influenced by the 
assumptions of the models, but also by the specifics of the observa- 
tions. The spatial coverage of the kinematics is one example. The 
two-dimensional coverage is an improvement over a few slits of- 
ten used in other studies. Similarly, increasing the radial extent of 
the data could change the results. Another issue, associated with 
the modelling techniques, is the ability of the three-integral models 
to recover the true distribution function of the galaxy. This is very 
important for the investigation of the internal dynamics, since the 
recovered orbital distribution must represent the observed galaxy if 
we want to learn about the galaxy's evolutionary history. In this pa- 
per we present tests designed to probe these issues and, in general, 
to determine the robustness of our three-integral method. 

This paper is organised as follows. Section|2|summarises the 
SAURON spectroscopy and the photometric ground- and space- 
based data. The analysis of the velocity maps, used to quantify the 
presence and influence of possible non-axisymmetric motions as 
well as a brief discussion on bars in NGC 2974, is presented in Sec- 
tion|3| The three-integral dynamical models for the stellar motions 
are discussed in Section|4| Section|5]is devoted to tests of the three- 
integral method involving the determination of the model parame- 
ters (T, j), influence of the radial extent of the data and the recovery 



Table 1. Properties of NGC 2974. 



Parameter Value 

Morphological type E4 

Ms [mag] -20.46 

effective B-V [mag] 1.00 

PA [degrees] 42 

Distance Modulus [mag] 31.66 

Distance scale [pc/arcsec] 104.13 



Notes - Listed properties are taken from the Lyon/Meudon Extragalactic 
Database (LEDA). Distance modulus is from .Tonrv et al. 1 200U 



of the DF. The modelling of the emission-line gas kinematics and 
comparison with the results of the stellar dynamical modelling is 
presented in Section|6| Section0concludes. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of NGC 2974 used in this work consist of 
ground- and space-based imaging, and ground-based integral-field 
spectroscopy. The imaging data were presented in EGF03 and 
the absorption- l ine ki nematics of the SAURON observations in 
'Emsellem et al.' ('2004, hereafter E04) as part of the SAURON sur- 
vey ide Zeeuw et al. 2002). In this study we also use the SAURON 
emission-line kinematics of NGC 2974. 



2.1 SAURON spectroscopy 

NGC 2974 was observed with the integral-field spectrograph 
SAURON mounted on the 4.2-m William Herschel Telescope 
(WHT) in IVIarch 2001. The observations consisted of eight ex- 
posures divided equally between two pointings, each covering the 
centre and one side of the galaxy. The individual exposures of both 
pointings were dithered to obtain a better estimate of detector sensi- 
tivity variations and avoid systematic errors. The instrumental char- 
acteristics of SAURON and a summary of the observations are pre- 
sented in Table|2| 

The SAUR ON data were reduced following the steps described 
in IB aeon et alj i2001.) using the dedicated software XSauron de- 
veloped at CRAL-Observatoire. The performed reduction steps in- 
cluded bias and dark subtraction, extraction of the spectra using 
a fitted mask model, wavelength calibration, low frequency flat- 
fielding, cosmic-ray removal, homogenisation of the spectral reso- 
lution over the field, sky subtraction and flux calibration. All eight 
exposures were merged into one data cube with a common wave- 
length range by combining the science and noise spectra using op- 
timal weights and (re)normalisation. In this process we resampled 
the data-cube to a common spatial scale (CK'S x O'.'S) with a resulting 
field of view of about 45" x 45" . The data cube was spatially binned 
to increase the signal-to-noise (S /N) ratio over the fie l d, usin g the 
Voronoi 2D binning algorithm of ICappellari & CopinI i2003 t). The 
targeted minimum 5'/A'^ was 60 per aperture, but most of the spec- 
tra have S/N ratio high (e.g. [S/N]r„ax ~ 420) and about half 
of the spatial elements remain un-binned. The final data cube of 
NGC 2974 and the detailed reduction procedure was presented in 
E04. 
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Figure 1. Distribution of H/3 and [OIII] emission lines and gas kinematics observed by SAURON. Gas intensities are in mag arcsec"^ witli arbitrary zero 
points. Gas mean velocity V and velocity dispersion a are in km s~^. Overplotted contours are levels of the reconstructed total intensity from the full SAURON 
spectra. Arrow in the left corner of each plot marks North-East orientation of the maps. 



Table 2. The SAURON instrumental characteristics and exposure details of 
the observations of NGC 2974 obtained in March 2001 at the WHT. The 
sky apertures are pointed 1 '.9 away from the main field. 



Field of view 


33" X 41" 


Aperture size 


0'.'94 


Final spatial sampling 


0"8 


Spectral range 


4810 - 5300A 


Spectral sampling 


1.1 A pixel" 1 


Spectral resolution 


4.2A (FWHM) 


# of field lenses 


1431 


# of sky lenses 


146 


# of exposures 


8 


Exposure time per pointing 


1800 s 


Instrumental dispersion (a) 


108 kms-l 


Median seeing (FWHM) 


1'.'4 



2.2 Absorption-line kinematics 

The SAURON spectral range includes several important emis- 
sion lines: H/3, [Oin]AA4959,5007 and [NI]AA5198,5200 dou- 
blets. These lines have to be masked or removed from the spec- 
tra used for the extraction of the stellar kinematics. The method 
most suitable for this is the direct pixel-fitting method operat- 
ing in wavelength space, which allows easy masking of the emis- 
sion lines. We used the penalis ed pixel-fitting algorithm (pPXF) of 
ICappellari & Emsellemli2bo4) . following the prescriptions of E04. 
The line-of-sight velocity distribution (LOSVD) was parametrised 
by the Gauss -Hermite expansion ( van der Marel & Franx 1993; 
lGerharcil99'3l) . The 2D stellar kinematic maps of NGC 2974, show- 
ing the mean velocity (V), the velocity dispersion (cr), as well as 
higher order Gauss-Hermite moments hs and hi, were presented in 
E04, along with the kinematics of 47 other elliptical and lenticular 
galaxies. In this study we expand on the previously published kine- 
matics by including two more terms in the Gauss-Hermite expan- 
sion (/i5 and he) to make sure all useful information was extracted 
from the spectra and tighten the constraints on the dynamical mod- 
els. The extraction of additional kinematic terms was performed 
following the same procedure as in E04. The new extraction is con- 
sistent with the published kinematics (except that now the LOSVD 
is parameterised with 6 moments) and we do not present them here 
explicitly (but see Fig. ll2t . 

We estimated the errors in the kinematic measurements by 



means of Monte-Carlo simulations. The parameters of the LOSVD 
were extracted from a hundred realisations of the observed spec- 
trum. Each pixel of a Monte-Carlo spectrum was constructed 
adding a value randomly taken from a Gaussian distribution with 
the mean of the observed spectrum and standard deviation given 
by a robust-sigma estimate of the residual of the fit to the ob- 
served spectrum. All realisations provide a distribution of values 
from which la confidence levels were estimated. During the ex- 
traction of the kinematics for error estimates, we switched off the 
penalisation of the pPXF meth od in order to obtain the true ( unbi- 
ased) scatter of the values (see lCappellari & Emselleri]l2004l for a 
discussion). 



2.3 Distribution and liinematics of ionised gas 

NGC 2974 has p reviously been se arched for the existence of 
emission-line gas. iKim et all il988l) reports the detection of HI 
in a disc structure aligned with the optical isophotes. The total 
mass of HI is estimated to be 8 x 10^ Mp, ro tating in a disc with 
an inclination of i ~ 55°. Bus on et alJ il993l) detected Ha emis- 
sion distributed in a flat structure along the major axis. Assum- 
ing a disc geometry, the inferred inclination is ~ 59°, and the 
total mass of HII was estim ated to be ~ 3x IC* M0. Similar re- 
sults are found bv lPlana'et a l. 1 1998). Deep optical ground-based 
imaging studies suggested the existence of "arm-like" spiral struc- 
tures, visible in the filamentary distribution of ionised gas outside 
~ 5" iBregman et al. 1992; Buson et al. 1993). The recent high- 
resolution HST imaging in Ha-I-[NII] revealed the presence of a 
gaseous two-arm spiral in the inner ~ 200 pc, with a total mass of 
6.8 X 10'' Mq (EGF03). 

The strongest emission line in the SAURON spectra of 
NGC 2974 is the [OIII] doublet. There is also considerable emission 
in H/3 and some emission from the [Nl] lines. Measurement of the 
emission-line kinematics followed the extraction of the absorption- 
line kinematics. For each spectrum in the data-cube we performed 
three steps: 

(i) The pPXF method provided the model absorption spectrum 
that yielded the best fit to the spectral range with the emission lines 
([OIII], H/3 and [Nl]) excluded. 

(ii) We then subtracted the model absorption spectrum from the 
original observed spectrum. This resulted in a "pure emission-line" 
spectrum which was used to extract the gas kinematics. 
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Table 3. Summary of the ground- and space-based observations of NGC 
2974. The exposure times of the HSTAVFPC2 observations are averages of 
all frames used to produce the WFPC2 association images. 





JKT 


HSTAVFPC2 


Filter band 


/ 


F547M&F814W 


Exposure time (s) 


60 


700 & 250 


Field of view (arcsec) 


380 X 350 


32 X 32 


Pixel scale (arcsec) 


0.3106 


0.0455 


Date of observations 


16.04.1993 


16.04.1997 



(iii) Each emission line was approximated with a Gaussian. The fit 
was performed simultaneously to the three lines of [OllI] and H/3, 
not using the mostly negligible [Nl] doublet. 

This procedure assumes that the velocity and velocity disper- 
sion of the different emission lines are equal. Performing the si- 
multaneous fit to the lines while allowing them to be kinematically 
indepe ndent yields similar results (Sarzi et al. in preparation). Fol- 
lowing lOsterbrockl fl989) we assumed a 1:2.96 ratio for the com- 
ponents of the [OIII] doublet, while leaving the intensity of the 
[OIII] and H/3 lines independent. The flux maps of [OIII] and H/3 
lines as well as the maps of the [OIII] emission-line mean velocity 
and velocity dispersion are presented in Fig. Q The instrumental 
broadening of 108 kms~^ was subtracted in quadrature from the 
emission-line velocity dispersion map presented here and also used 
in Section|6| 

Both the H/3 and [OIII] emission lines are present over the 
whole extent of the maps on Fig.Q but their intensity drops off 
approximately exponentially with distance from the centre. The 
[OIII] emission is stronger over the entire SAURON field with the 
[OIII] to H/3 line-ratio being ~1.7. The shape of the distributions 
are very similar, although H/3 follows the stellar light isophotes 
more precisely. The [OIII] distribution shows departures from the 
stellar isophotes in two roughly symmetric regions, positioned east 
and west of the centre on Fig. Q The nature of these dips in the 
[OIII] flux are discussed in Section l33l 

2.4 Ground- and space-based imaging 

In this study we used the existing ground- and space-based images 
of NGC 2974. The already reduced wide-field ground-based /-band 
image of NGC 2974 was taken from tGoudfrooii et al. ( J994), ob- 
tained at the 1.0-m Jacobus Kapteyn Telescope (JKT). We also re- 
trieved the Wide Field and Planetary Camera 2 (WFPC2) associa- 
tion images of NGC 2974 from the Hubble Space Telescope (HST) 
archive (Program ID 6822, PI Goudfrooij). The details of all imag- 
ing observations are presented in Table|3| 

A major complication in the derivation of the surface bright- 
ness model needed for the dynamical modelling is the existence 
of dust, clearly visible on the high resolution images. We consid- 
ered two possible approaches: masking the patchy dust areas and 
excluding it from the calculation of the model, or constructing a 
dust-corrected image. We decided to adopt the latter approach to 
determine the stellar surface brightness. We derived the correction 
of dust absorption using the F547M and F81 4W WFPC2 images, 
following the steps listed in ICappellari et"ai] |2002). The process 
consists of construction of a colour excess map E(V-I), from a cal- 
ibrated V-I colour image. The colour excess map is used to correct 
the pixels above a given E(V-I) threshold using the standard Galac- 
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Figure 2. V-I colour versus the elliptical radius of every pixel in the inner 
15" of the WFPC2/PC1 images of NGC 2974. The straight hne presents the 
best fit to the points obtained by minimising the absolute deviation. Notice 
the excess of red pixels between 0'.'3 and 2'.' 5 caused by dust. 



tic extinction curve. We assumed that the dust is a screen in front 
of the galaxy and that dust-affected pixels have the same intrinsic 
colour as the surrounding unaffected pixels. Figure|2|shows the cal- 
ibrated V-I colour of pixels in the inner part of the PC images. The 
best fit to the colours was obtained by minimising the absolute de- 
viation of the pixel values. This fit, represented by a line in Fig.|2| 
was used to calculate the colour excess by subtracting the measured 
colour from the fit. The resulting E(V-I) image is shown in the sec- 
ond panel of Fig. |3| The other panels on the same figure present 
the inner parts of the F814W PC image before and after the cor- 
rection of dust absorption. The colour excess image highlights the 
dust structure visible also on Fig. 3 of EGF03 and suggests a non- 
uniform distribution of dust in the central region of NGC 2974. 



3 QUANTITATIVE ANALYSIS OF VELOCITY MAPS 

Two-dimensional kinematic maps offer a large amount of informa- 
tion and are often superior to a few long-slit velocity profiles. The 
two-dimensional nature of these data motivates us to quantify the 
topology and structure of these kinematic maps, just as is com- 
monly done for simple imaging. We have developed a new tech- 
nique to deal with kinematic maps based on the Fourier expansion, 
and, due t o its si milari ty to the photometry, we named it kineme- 
try ilCoDinetal][200lh . This method is a generalis ation of the 
approach devel oped for tw o-dimensional radio data iFranx et'al] 
1994; Schoenm akers et all 1997: Wong et al. 20o3). The aim of the 
method is to extract general properties from the kinematic maps of 
spheroidal systems (early-type galaxies) without assuming a spe- 
cific intrinsic geometry (e.g. thin disc) for the distribution of stars. 
This changes the interpretation and the approach to the terms of the 
harmonic expansion from the case of cold neutral hydrogen or CO 
discussed in the above-mentioned papers. In this section we briefly 
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Figure 3. Dust correction on the F814WWFPC2 image of NGC 2974. From left to right: observed F814W image; colour excess E(V-I) obtained as described in 
the text; dust corrected F814W image. The arrow points to the north and associated dash to the east. All images were constructed using histogram equalisation. 
Lighter shades represent brighter regions. The E(V-I) map is stretched between -0.1 (dark) and 0.2 (bright) magnitudes. 



present the method and apply it to the stellar and gaseous velocity 
maps (Krajnovic et al. (in prep.)). 

3.1 Harmonic Expansion 

The kinemetry method consists of the straightforward Fourier ex- 
pansion of the line-of-sight kinematic property K{r, 9) in polar co- 
ordinates: 

JV 

K{r,e) = ao(r) + ^ c„(r-) cos[n(e - 0„(r))]. (1) 

n = l 

The expansion is done on a set of concentric circular rings (al- 
though other choices are possible), and its main advantage is lin- 
earity at constant r. The expansion is possible for all moments of 
the LOSVD, but in this paper we restrict ourselves to the mean ve- 
locity maps. 

The kinematic moments (moments of LOSVD) of triaxial 
galaxies in a stationary configuration have different parity, e.g., 
mean velocity is odd, while the second moment (u^), is even. The 
parity of a moment generates certain symmetries of the kinemat- 
ics maps. More generally, the maps of odd moments are point — 
anti — symmetric, or: 

V{r,e + n) ^ -V{r,e). (2) 

If axisymmetry is assumed, in addition to the previous relation, 
maps are mirror — anti — symmetric, or: 

V{r,7v-e)^-V{r,e), (3) 

These symmetry conditions translate into the requirement on the 
harmonic expansion (eq.0 that for point-anti-symmetric maps the 
even coefficients in the expansion are equal to zero, while in the 
case of mirror-anti-symmetry, additionally, the odd phase angles 
have a constant value, equal to the photometric position angle (PA) 
of the galaxy in the case of a true axisymmetric galaxy. This means 
that to reconstruct the mean velocity map of a stationary triaxial 
galaxy, it is sufficient to use only odd terms in the expansion. 

These properties of the velocity maps enable certain natu- 
ral filtering (point-(anti)-symmetric — eq. J2j, and mirror-(anti)- 
symmetric — eq. js}) using the harmonic expansion with coeffi- 
cients set to zero or phase angles fixed at certain values. For (visual) 
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Figure 4. Kinemetric expansion of the stellar velocity field as a function of 
radius. From top to bottom: (a) first amplitude coefficient in the harmonic 
expansion ci; (b) ratios of amplitude coefficients presented by: triangles 
C2/C1, diamonds C3/C1, squares C4/C1 and plusses cg/ci; (c) important po- 
sition angles: phase angle ij>i representing the kinematic position angle (di- 
amonds), photometric position angle as measured on WFPC2/PC1 F814W 
image (asterisks) and adopted value for photometric position angle mea- 
sured from the integrated SAURON flux image (dashed straight line); (d) 
difference between the first and third phase from the kinemetric expansion. 



comparisons of the data with the results of axisymmetric modelling 
it is useful to apply the axisymmetric filtering to the data, as we will 
see below (Section|6). 
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Figure 5. Same as Fig|3 but for the gas velocity map. 



3.2 Kinemetric analysis of velocity maps 

We wish to know the intrinsic shape of NGC 2974 and, in par- 
ticular, whether it is consistent with axisymmetry, which would 
permit the construction of three-integral axisymmetric dynamical 
models of the galaxy. In order to obtain the necessary information 
we applied the kinemetric expansion to the observed stellar veloc- 
ity map. If NGC 2974 is an axisymmetric galaxy, the kinemetric 
terms should have odd parity (even terms should be zero) and the 
kinematic position angle should be constant and equal to the PA. 

The amplitude and phases of the first five terms in the expan- 
sion are presented in Fig. |4| The first panel presents the dominant 
term in the expansion, ci , which gives the general shape and ampli- 
tude of the stellar velocity map. The correction to this term is given 
by the next significant term, C3, which is already much smaller 
than ci and is presented in the second panel as a fraction of ci. 
Even terms in the expansion, C2 and C4, are also plotted on the 
same panel, and are much smaller (~ 1% of ci)^. For compari- 
son, the fifth term in the expansion, C5, is also plotted and is larger 
than both C2 and C4. The SAURON pixel size is ~ 1" and mea- 
surement of higher-order terms at radii smaller than 2" cannot be 
trusted. Clearly, the velocity map in NGC 2974 can be represented 
by the first two odd terms in the expansion. Neglecting all higher 
terms results at most in a few percent error. 

The lower two panels in Fig.|4|present the phases of the dom- 
inant terms. The cj>\ phase is defined as the kinematic angle of the 
velocity map, here measured east of north. This angle is compared 
with measurements of two important angles: (i) the PA measured 
on the WFPC2/PC F814W dust-corrected image using the I RAF 
ellipse fitting task ellipse and (ii) the PA measured on the recon- 
structed SAURON flux image obtained by integrating the spectra in 
each bin. The agreement between the different angles measured on 
the SAURON observations is excellent, with slight departures in the 



^ The zeroth term, ao , gives the systemic velocity of the galaxy and is not 
important for this analysis. 
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Figure 6. Equivalent widths of [OIII] and H/3 emission-lines. These plots, 
constructed by dividing the emission-line flux by the symmetrised stellar 
continuum, show the relative distribution of emission-lines to the stellar 
continuum. The dashed ellipse represents the possible OLR ring of the inner 
bar (see text for details). Overplotted contours are levels of the reconstructed 
total intensity from the full SAURON spectra. The arrow in the left comer of 
each plot marks North-East orientation of the maps. 



inner 3". The PA measured on the high resolution WFPC2 image 
suggests a small photometric twist in the inner 10" of « 3°, which 
can be also seen in Fig. llOl 

The phase angle (j)z is the phase of the third term in the kine- 
metric extraction. It is easy to show, if the galaxy is axisymmetric 
(requiring in eq.Q K(r,^) = for 6 — 4>i + 7r/2), and the higher 
terms can be neglected, that the phases (jii and ^3 satisfy the rela- 
tion: 



nil 



(4) 



where n G Z. The last panel in Fig.|4]shows this phase difference. 
The condition given by eq.(|4j is satisfied along the entire investi- 
gated range with a small deviation in the inner 3". Summarising all 
the above evidence, we conclude that the observed stellar kinemat- 
ics in NGC 2974 is consistent with axisymmetry. 

We repeated the kinemetric analysis^ on the emission-line gas 
velocity maps and present the results in Fig.|5] While the amplitude 
coefficients, Ci, are similar to the stellar coefficients (small values 
of all terms higher than C3), the behaviour of the phase angles is 
quite different. 

The last panel of Fig.|5|is perhaps the best diagnostic tool. The 
dashed line presents the required value for the difference between 
the phase angles, (f)i — tp^, assuming axisymmetry. Deviations are 
present in the inner 4" and, although much smaller, between 9" 
and 11". These deviations indicate departures from axisymmetry, 
which are strongest in the central few arcsecs. 

3.3 Signature of bars in NGC 2974 

In the previous section, we quantified the signatures of non- 
axisymmetry on the gas velocity maps, especially strong in the 
inner 4". Inside this radius, EGF03 discovered a two-armed spi- 
ral and explained it by a weak bar with the corotation resonance 
(CR) at 4'.'9 and outer-Lindblad resonance (OLR) at 8'.'5. These 
scales are consistent with the observed departures from axisym- 
metry in the SAURON kinematic maps. An additional confirmation 

^ A similar analysis app roach for a gas disc would be using the 
*Scho enmakers etalJiigg T'ihar monic analysis on a tilted-ring model of the 
gas di sc, interpreting the results within epicycle theory (see also Wong et al] 
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Figure 7. Diagram of resonances in NGC 2974, derived from the po- 
tential of the best-fitting stellar dynamical model. The upper horizontal 
line is the pattern speed of the nuclear (secondary) bar from the EGF03 
study (Qp = 700 kms^^ kpc~^). The lower horizontal line shows 
the position of the inferred pattern speed of the lai'ge-scale primary bar 
(Op = 185 km s~^fcpc~^). The vertical lines show the assumed radial 
positions of important resonances of the primary bar. From left to right: 
ILR (CR of the secondary bar), UHR, CR and OLR. 



comes from the equivalent width maps of emission-lines (Fig.|6l. 
These maps show the distribution of the emission-lines relative to 
the stellar continuum. Especially noticeable is the ring structure in 
the [OIII] equivalent width distribution. For comparison, we over- 
plotted an ellipse with semi-major axis length of 8'.'5 (the radius of 
OLR of the EGF03 bar). The orientation (~ 25° away from the ma- 
jor axis of the galaxy) and the size of the ring ('^ 8'.' 5) are in good 
agreement with the expected characteristics of the bar discovered 
by EGF03. 

The equivalent width maps show that [OIII] emission-line is 
also influenced outside the OLR. There is an increase in the value 
of the equivalent width at larger radii as well as a filamentary (spi- 
ral) structure connecting the ring and this large scale region. Also, 
as seen before on Fig.Q the distribution of the [OIII] emission-line 
intensity exhibits an elongated structure in the central 5". Simi- 
larly, around 8" — 10" [OIII] is also elongated, but this time ap- 
proximately perpendicular to the first elongated structure. Follow- 
ing this, the [OIII] map has a plateau between 12" and 15". The 
end of the plateau is followed by a dip in the distribution with a 
possible turn up at radii larger than 20". 

Although the central structure and kinematics of the [OIII] 
distribution are influenced by the inner bar, the large-scale struc- 
ture (beyond ~ 10") is not likely to be influenced by this weak 
inner bar. On a more speculative basis, we ca n infer the exis- 
tence of a large-scale bar. Studies of double bars I'Erwin & Soarkd 
L2OO2; Laine et al. 2002; Erwin & Soarke 2003) suggest a size ra- 
tio of about 5 to 10 between the primary and secondary bars; this 
would set the size of an hypothetical large-scale (primary) bar in 
NGC 2974 to between about 12'.' 5 and 25". 

A more precise, although still approximate, estimate of the 



properties of the primary bar can be obtained from the resonance 
curves of NGC 2974. Using the potential of the best-fitting stellar 
dynamical model (see Section|4} we constructed the resonance dia- 
gram presented in Fig.Q We calculated profiles of f2, — k/2, 
Q + k/2 and Q — k/4, where Q is the angular velocity, V/r, 
and K is epicyclic frequency, defined as = — h 4f2^. As- 
suming that the resonances of the primary and the secondary bar 
are coupled to m inimise th e chaos produced arou nd the resonances 
( iFriedli & 1^^1^161.1993; iPfenniger & Normaijfl99a) . where the 
inner Lindblad resonance (ILR) of the primary bar is the CR of the 
secondary bar, we can estimate the main resonances of the primary 
bar (indicated by the vertical lines on Fig.0. The pattern speed 
of the primary component is then ~ 185 kms~^ kpc~^, ILR is at 
4'.' 9, ultra harmonic resonance (UHR) at ^ 12'.' 5, CR at « 18'.' 6 
and OLR at ~ 27'.' 7. The size of the primary bar can be taken to be 
80% of its CR, so about 13 to 14"(~ 1.4 kpc). 

The presented analysis is strictly valid only for an axisym- 
metric potential with a weak bar perturbation, hence the above- 
mentioned size estimates are only approximate, but indicative. We 
see several features in the [OIII] distribution and velocity maps that 
support the assumption of a large-scale primary bar: the dip in the 
gas velocity map around 12", the plateau between 12"and 15"in 
the [OIII] distribution as well as the dip in the [OIII] distribution 
around 18". The last one corresponds to the position of the CR, 
which is a chaotic region devoid of gas, consistent with the ob- 
served lower flux in that region. 



3.4 Case for axisymmetry in NGC 2974 

The alignment of the gaseous and stellar component, previously 
detected and also confirmed in this study, suggests that NGC 2974 
is an axisymmetric galaxy. However, the gaseous component shows 
signatures of non-axisymmetric perturbations. The contribution of 
the non-axisymmetric motion, AV/ci, to the total velocity field 
can be quantified from the phase difference (/)3 — . If the condition 
in eq. ^4} is not satisfied then: 

= —sm3{(j)i - (j)3), (5) 

Cl Ci 

which is presented in Fig.|8|for both stellar and gaseous velocity 
maps. At ~ 3", AV/ci for the emission-line gas is ~ 0.1 and at 
~ 10" it is « 0.04, confirming that the non-axisymmetric contribu- 
tion is significant in the centre of the emission-line velocity map. Its 
influence on the stellar velocity is not significant over the SAURON 
field. Emission-line gas is a more responsive medium and unlike 
the stars, due to the viscosity of the gas particles, shows evidence 
of weak non-axisymmetric perturbations. It is possible that other 
early-type galaxies with disc-like components harbour such weak 
and hidden bar systems. 

Summarising, the stellar velocity map is point- and mirror- 
anti-symmetric, supporting an axisymmetric shape for NGC 2974. 
On the other hand, the gaseous velocity map shows strong devia- 
tions from mirror-anti-symmetry in the centre, and the distribution 
and equivalent width of [OIII] emission lines supports the weak in- 
ner bar found by EGF03 and suggests the existence of a weak large- 
scale bar. However, since the bar perturbations on the axisymmetric 
potential are presumed to be weak and appear do not influence the 
stellar kinematics, we ignore them in the remainder of the paper, 
and describe NGC 2974 with an axisymmetric potential. 
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Figure 8. Contribution of the non-axisymmetric motion to tlie stellar and 
gaseous velocity fields, as a fraction of the dominant term in the kinemetric 
expansion, c\ . Diamonds represent the emission-line gas velocity contribu- 
tion, while triangles denote the stellar velocity contribution. 



Table 4. The MGE parameters of the circular PSF of HSTAVFPC2/F814W 
filter. 



Table 5. The parameters of the MGE model of the deconvolved I-band 
WFPC2 surface brightness of NGC 2974. Columns present (from left to 
light): number of the two-dimensional Gaussian, central intensity of the 
Gaussian, width (sigma) of the Gaussian, axial ratio of the Gaussian, total 
intensity contained in the Gaussian. 
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4 STELLAR DYNAMICAL MODELLING 

In order to investigate the orbital structure of NGC 2974 we con- 
struct fully general axisymmetric models of the galaxy's stellar 
component. The three-integral models prese nted here are based on 
Schwarzschild's orbit superposition method fechwarzschilt^' 19791 
M^SA, furth er develo ped by Rix et al. ( 1997), van der Marel et al. 
i 19981) and ICretton et al. 1 1999), and adapted for more general 
surface-brightn ess distributions by Cappellari et aL (2002, here- 
af ter C02') and [^rolm^t^ hereafter V02), similarly as 

in ICretton & van denBoschH l999'). The three-integral modelling 
technique is widely used for constructing dynamical models of ax- 
isymmetric galaxies, and has been thoroughly described in the lit- 
erat ure by the aforementioned authors as we ll as by other groups 
(e.g. lOebhardt et alj2003tlvalluri et alJ2o'ol . It is most commonly 
used to determine the masses of the central black holes in nearby 
galaxies and investigate the internal orbital structure of the galax- 
ies. The SAURON observations of NGC 2974 do not have the neces- 
sary resolution to probe the sphere of influence of the central black 
hole and we therefore restrict ourselves to the determination of the 
mass-to-light ratio T and the inclination i of the galaxy, as well as 
the internal orbital structure. 

4.1 The Multi-Gaussian expansion mass model 

The starting point of the stellar dynamical modelling is the deter- 
mination of the gravitational potential of the galaxy. The potential 
can be obtained by solving the Poisson equation for a given den- 
sity distribution which can be derived by deprojecting the obser- 
vations of the 2D stellar surface density. I n this work we used th e 
multi-Gaussian expansion (MGE) method jEmsellem et alj[l994h . 
following the approach of C02 and V02. 

In order to get the MGE model, we simultaneously fitted the 
ground-based I-band image and the dust-corrected PC part of the 



W FPC2/F814W ima ge using the method and software developed 
bv ICappellaril i2002h . The dust correction (Section l2.4t success- 
fully removed the dust contamination from the high-resolution im- 
age of the nucleus, but the large-scale image was badly polluted by 
several stars, with a particularly bright one almost on the galaxy's 
major axis. We masked all stars inside the model area to exclude 
them from the fit. The ground-based image, used to constrain the fit 
outside 25", was scaled to the WFPC2/PC1 image. We computed 
the PSF of the F814W PCI image at the position of the nucleus 
of NGC 2974, using the TinyTim software (Krist & Hook 2001), 
and parametrised it by fitting a circular MGE model with constant 
position angle as in C02. Table |4| presents the relative weights G-, 
(normalised such that their sum is equal to unity) and the corre- 
sponding dispersions a of the four Gaussians. Table |5| gives the 
parameters of the MGE model an alytically deconvo lved from the 
PSF. Following the prescription of lCappellaril i2003) . we increased 
the minimum axial ratio of the Gaussians, qj, until the signif- 
icantly changed, in order to make as large as possible the range 
of allowed inclinations by the MGE model. The upper limit to the 
q-j was also constrained such that the MGE model is as close as 
possible to a density stratified on similar ellipsoids. Although the 
d eprojection o f an axisymmetric density distribution is non-unique 
( Rvbicki 1987), our 'regularisation' on the MGE model produce re- 
alistic intrinsic densities, while preventing sharp variations, unless 
they are required to fit the surface brightness. We verified that the 
MGE model used in this study is consistent with the MGE model 
presented in EGF03. 

The comparison between the MGE model and the photome- 
try along different angular sectors is shown in Fig.|9| The profiles 
are reproduced to within 3%, and the RMS error is about 2%. The 
increase of the relative error (right-hand panel in Fig. |9) at larger 
radii is caused by the light pollution of the bright star on the ma- 
jor axis of the galaxy and can be better understood by looking at 
the comparison of the convolved model and the actual observa- 
tion. Fig. llOl presents both ground- and space-based images and the 
MGE model. On the WFPC2/PC image there is a slight deviation 
from the model (of constant PA) about 10" from the nucleus. The 
deviating structure is point-symmetric and r eminiscent of a spiral 
perturbation , as su ggested by other studies iBresman et aljfT993 
Buson et al] | 19931 : EGF03). Except for this slight departure from 
the model the galaxy surface brightness is well represented by the 
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Figure 9. Left Panels: Comparison between the combined WFPC2/F814W 
and ground-based photometry of NGC 2974 (open squares) and the MGE 
model (solid line) in seven angular sectors as function of radius. The indi- 
vidual convolved Gaussians are also shown. Right panels: radial variation 
of the relative error along the profiles. 

MGE model and we use it to calculate the representative gravita- 
tional potential. 

4.2 Construction of three-integral models 

Briefly, Schwarzschild's method can be divided in four steps. In 
the first step, one derives the stellar potential assuming the shape 
of the potential (axisymmetry) and stellar mass-to-light ratio, T 
(free parameter), through deprojection of a parametrisation of the 
surface density (in this case MGE parametrisation which can be de- 
projected analytically). The second step involves the construction 
of a representative orbit library by integrating the orbits in the de- 
rived potential. Each orbit is specified by three integrals of motion 
(E, Lz, I3), where E is the energy, Lz the component of the angu- 
lar momentum along the z symmetry axis, and Is is a non-classical 
integral, which is not known analytically. The integral space is con- 
structed on a grid that includes >99% of the total luminous mass of 
the galaxy. The next step consists of projecting the orbits onto the 
space of observables (x', y', Vtos), where (x' , y') are in the plane 
of the sky and vios is the line-of-sight velocity given by the obser- 
vations. In our implementation, this is done taking into account the 
PSF convolution and aperture binning. The final step of the method 
is to determine the set of weights for each orbit that, when added 
together, best corresponds to the observed kinematics in the given 
spatial bin as well as reproducing the stellar density. In our imple- 
mentation of the method, this best-fitting set is found by solving a 
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Figure 10. Contour maps of the ground-based I-band and dust-corrected 
WFPC2/F814W images. The brightest star on the ground-based image as 
well as four additional stars covered by the model were masked out and 
excluded from the fit. Superposed on the two plots are the contours of the 
MGE surface brightness model, convolved with the WFPC2 PSF. 

non-negative least-squares (NNLS) problem using the routine writ- 
ten by Lawson & Hanson ( 1974). 

The software implementation used here is similar to that used 
in the V02 and C02 studies, but it has evolved substantially since. 
The improvements are described in detail in Cappellari et al. (in 
prep.). We verified that the results from the new code are the same 
as from the old one if identical settings are adopted. An app lication 
to the elliptical galaxy NGC 4473 was presented in Cappell ari et alJ 
Here we give a quick overview of the changes with respect 
to the description in C02: 

1. The method requires that the orbits sample a three- 
dimensional space of integrals of motion, the energy E, Lz and J3. 
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Figure 11. A grid of inclination angle i versus M/L ratio T. Contours 
present constant Ax^, measuring the goodness of fit of the dynamical mod- 
els. Every symbol corresponds to a three-integral axisymmetric model with 
given inclination angle and T ratio. The best-fitting model for each incli- 
nations are presented with grey symbols. The first three contours show the 
formal 68.3%, 95.4% and 99.73% (thick grey contour) confidence regions 
for two degrees of freedom. The best-fitting model on the grid has an incli- 
nation of 65°and T = 4.5. 

In the new scheme fsee lCretton et alll999l for details of the previ- 
ous approach), at each E, we construct a polar grid of initial start- 
ing positions on the meridional plane (linear in angle and in radius), 
going from i? = z = to the curve defined by the thin tube orbits 
(to avoid duplication) which is well approximated by the equation 
E? + = R^{E), where Rc{E) is the radius of the circular orbit 
at energy E. The orbits are released with vr = Vz = and 0. 
In this way we sample the observable sp ace by uniformly distrib ut- 
ing the position of the orbital cusps (see lCappellari et all2004^ on 
the sky plane. 

2. Improved treatment of seeing effects and instrumental point 
spread function (PSF) by a Monte Carlo method. The PSF can be 
considered as the probability that an observed photon arriving at 
a detector will be displaced from its original position by a given 
amount (specified by the PSF characteristics). The projected or- 
bital points (results of the orbit integration and projection onto the 
sky-plane) are stored in the apertures in which they landed after ap- 
plying a random displacement taken from the Gaussian probability 
distribution defined by the PSF. 

3. Generalisation of the projection of the orbits into the space of 
the observables. The bins of the optimal Voronoi binning of the 
two-dimensional integral-field data have non-rectangular shapes. 
The orbital observables now can be stored on apertures of any shape 
that can be represented by polygons. 

4.3 Stellar dynamics - modelling results and discussion 

Our stellar kinematic maps of NGC 2974 consist of 708 Voronoi 
bins. Each bin contributes with 6 kinematic observables to which 
we also add the intrinsic and projected mass density observables. 



resulting in a grand total of 5664 observables. The largest orbit li- 
brary that was computationally possible for the given number of 
observables, consists of 2 x 41 x 10 x 10 = 8200 orbits (for each 
of the 41 different E we construct a polar grid of starting points 
sampling 10 angles and 10 radii). With this choice of orbit library, 
the number of observables is smaller than the number of orbits, and 
the NNLS fit will not have a formally unique solution. Moreover the 
recovery of the orbital weights for the orbits from the observations 
is an inverse problem, and as such is intrinsically ill-conditioned. 
For these reasons a direct solution of the problem generally con- 
sists exclusively of sharp isolated peaks. It is unlikely for the DF 
of real galaxies to be very jagged, since (violent) relaxation pro- 
cesses tend to smooth the DF. Moreover observational constraints 
on the smoothness of the DF, at least for the bulk of the stars in a 
galaxy, come from the smoothness of the observed surface bright- 
ness, down to the smallest spatial scales sampled by HST. 

A standard mathematical approach to solve inverse problems 
is by regularising (e.g. Press eta ij ii99a Chapter 18). This has 
been generally applied by all groups inv olved in this modelling 
approach (e.g. iRix et alJll997t Ib ebhardt e tai]|2003l : IValluri et al] 
l2004ICretton' ^^mseUM32004l) . Regularisation inevitably biases 
the solution, by forcing most orbits to have a non-zero weight. The 
key here is to apply the right amount of regularisation. Previous 
tests with the Schwarzschild code suggested a value of the regu- 
larisation parameter A = 4 (see Cretton et al. 1999; V02). After 
initial testing we also adopted this value. For more details on regu- 
larisation, see McDermid et al. (in prep.). 

The models are axisymmetric by construction. In order to 
avoid possible systematic effects we additionally symmetrise the 
stellar kinematics as usually done in other studies. The sym- 
metrisation uses the mirror-(anti)-symmetry of the kinematic 
fields, such that kinematic values from four symmetric positions 
({x,y), {x, —y); {—X, y), {—x, —y)) were averaged. However, the 
Voronoi bins have irregular shapes and they are not equally dis- 
tributed with respect to the symmetry axes of the galaxy (minor 
and major axes). In practice, we average the four symmetric points, 
and, if for a given bin there are no bins on the symmetric positions, 
we interpolate the values on those positions and then average them. 
As the number of data points is not decreased in this way, the errors 
were left unchanged. 

Finally, we fitted axisymmetric dynamical models to sym- 
metrised observations of velocity, velocity dispersion and Gauss- 
Hermite moments (/is — ha) while varying the mass-to-light ra- 
tio T and the inclination i. Figure fTTI shows a grid of our mod- 
els with overplotted Ax^ contours. The best-fit parameters are 
T = 4.5 ± 0.1 (in the I band) and i = 65° ± 2.5°. The data - 
model comparison for these values is given in Fig. ll2l (svmmetrised 
data are shown in the first row and the best-fitting model in the 
third row). In the same figure we present the best-fitting models for 
the given inclinations. The formal statistical analysis firmly rules 
out (with 3cr confidence) all inclinations outside i — 65° ± 2.5°. 
The differences between the models are only marginally visible 
on Fig. 1121 However, comparing the models with the top row of 
symmetrised data, it is noticeable that the velocity dispersion is 
less well fitted with increasing inclination. On the other hand, the 
fit to /t3 improves with higher inclinations. Higher-order moments 
change similarly and the final is the result of this combined ef- 
fect. Surprisingly, the difference between the best-fit model and the 
data are bigger than the difference between the other models and 
the best-fit model, provoking a question whether the determination 
of the inclination is robust. 

Recently, EGF03 used two-integral f{E, Lz) dynamical mod- 
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Figure 12. Comparison of the symmetrised observation of NGC 2974 (first row) and four orbital superposition models with best mass-to-light ratio T, for 
a given inclination i. From second row to bottom: (i, T) = (55, 4.4), (65, 4.5), (75, 4.7), (89, 4.8). From left to right each panel presents: mean velocity V, 
velocity dispersion cr, and Gauss-Hermite moments: hs, hi, hs and feg. Isophotal contours of total light are shown with solid lines. 



els and found a best fit for an inclination of 60° but also stated that 
models with 58° < i < 65° fit equally well. Our three-integral 
models have a larger freedom in fitting the observations and it is 
unknown whether these models can uniquely constrain the incli- 
nation. Our obtained inclination is close to the previous measure- 
ments in the literature (references in SectionsQand|2j, and also to 
the inclination measured from gas and dust. NGC 2974 is perhaps 
a special case (in terms of its intrinsic structure and geometry) for 
which three-integral models are able to give a stronger constraint 
on the inclination. We return to this issue in Section l5!2l 

The integral space of the best-fitting model, (i.e. the space 
defined by the isolating integrals of motion (E, Lz,Iz) that define 
the orbits and the DF / = f{E,Lz,l3)), is shown in Fig. 1131 
Each panel presents mass assigned to orbits of constant energy, 
parameterised by the radius of the circular orbit. This orbit also has 
the maximum angular momentum and circular orbits of negative 
and positive angular momentum are in the bottom left and right 
corner of each plot, while the low angular momentum orbits are 
close to the symmetry (j/) axis. An interesting feature dominates 
the panels with radii larger then 5" 11. A high fraction of mass is 
assigned to orbits with high angular momentum. This indicates 
that the bulk of the stars between these radii rotate with high 



angular momentum. A possible physic al interpretation is that a 
large fraction of the stars orbit in a disc. lCinzano & van der Marell 
( 1994) argued that NGC 2974 has an embedded stellar disc and 
in their two-integral Jeans models they were not able to fit the 
stellar kinematics without introducing a disc with ~ 7% of the 
total galaxy light. On the other hand, the two-integral models of 
EGF03 did not need to invoke an additional stellar disc component 
to fit data consisting of their TIGE R data and the three long slits 
of ICinzano & van der Marel il994l) . The integral space presented 
here suggests that the three-integral models need orbits with high 
angular momentum, but the selected orbits also have different 
values of Li, and therefore do not represent a very thin stellar disc 
as assumed by Cinzano & van der Marel 1 1994), but a somewhat 
flattened distribution of stars like in a normal SO. The relative 
light contribution of the high-angular momentum orbits is ~ 10%, 
which corresponds to a total stellar mass of 1.5 X 10^° M0 
assuming the best-fit model inclination and T. 



1 2 Davor Krajnovic et al. 





":^m^.. '^.^^H^. "^^M^ 




Figure 13. Integral space of the best-titting orbit superposition model for NGC 2974. Each panel plots the meridional plane (R,z) with the starting positions 
of orbits (dots) for the given energy. Orbital starting positions correspond to the position of cusps (vji = =0). Overplotted is the fraction of mass assigned 
to orbits at constant energy, labelled by the radius of the circular orbit in arcseconds (printed at the upper left corner of each panel). The radius of the circular 
orbit is the size of the horizontal axis (measured from the centre). Negative values are obtained by sign reversal, since orbits can be prograde and retrograde. 
Orbits with high angular momentum are found in the right and left comer respectively on the plots. We show only the radii constrained by the data. The last 
three panels have the size of the SAURON field overplotted for comparison (red rectangle). The area inside the rectangle is constrained by the kinematic data. 
The colour bar on the right represents the relative colour coding of the orbital mass weights (7), where 1 represent the largest mass weight assigned to an orbit 
at the given energy. 



5 TESTS OF SCHWARZSCHILD'S 
ORBIT-SUPERPOSITION MODELS 

In the previous section we used three-integral models to recover 
the inclination, mass-to-light ratio and the internal structure of 
NGC 2974. Surprisingly, we find that the inclination is tightly con- 
strained by contours (Fig. ll H . although the difference between 
the models (Fig. I12> are smaller than the difference between for- 
mally the best-fitting model and the data. In this section, we wish 
to test the robustness of those results as well as the general abil- 
ity of the three-integral models to recover the given parameters. 
For this purpose we constructed an axisymmetric model mimick- 
ing NGC 2974 using two integrals of motion: the energy E, and 
the z-component of the angular momentum, ■ This two-integral 
galaxy model has the advantage of a known DF, / = f{E,Lz), ev- 
erywhere, which we want to compare with the results of the three- 
integral modelling. There are three issues we wish to test: 

1. Recovery of the input parameters of the two-integral model 
galaxy. This is a general test to show whether the three-integral 
method can recover the parameters used in construction of a test 
model. We wish to be consistent with the observations and consider 
only the recovery of the input mass-to-light ratio and inclination, 
especially in light of the results from the Section|4] 

2. Influence of the spatial coverage. The SAURON kinematic ob- 
servations of NGC 2974 roughly cover one effective radius. This is 
also the typical size of most kinematic observations of other early- 
type galaxies from the SAURON sample. Here we want to test the 
influence of the limited extent of the kinematic coverage on the re- 
covery of the orbital distribution. We do this by comparing the dif- 
ference between models using a limited and a. full spatial kinematic 
information provided by the two-integral galaxy model. 

3. The recovery of the input DF. We wish to test the abil- 
ity of our three-integral models to correctly recover the true 
input (two-in t egral) DF. Similar tests were also presented by 
iThomas et alJ i2004h for their implementation of Schwarzschild's 



method . 'Cretton et al.'l'l999|) and lVerolme & de Zeeuwl f2002^ and 
ICretton & Emsellemi 120041) described similar tests using two- 
integral Schwarzschild models. 



5.1 The input two-integral test model 



An f(E,L,) model 



model of NGC 2974 was constructed using the 
iHunte r & Oian 1 1993) contour integration method (hereafter HQ, 
see also Oian et al. 1995) . We used the mass model from Sec- 



tion 14. 1 1 parameterised by MGE. This approach follows in detail 
lEmselleme7ai]jl999h . and, due to the properties of Gaussian func- 
tions, simplifies the numerical calculations significantly. Given the 
density of the system, the HQ method gives the unique even part of 
the DF, /e = I [f{E, L,) + f{E, -L,)] (even in L^). The odd part 
can be calculated as a product of the even /e and a prescribed func- 
tion h = h{L/Lz). The magnitude of the function h is chosen to 
be smaller than unity, which ensures that the final DF, / — fs + fo, 
is physical (i.e. non-negative, provided that >> everywhere). 
In practice, the odd part is chosen to fit the observed kinematics 
(mean streaming) by flipping the direction of orbits with respect to 
the symmetry axis (photometric minor axis). 

The two-integral model of NGC 2974 was computed using 
i = 60° and T = 4.6. In order to construct a realistic DF we 
also included a black hole with mass Mbh ~ 2.5 x 10* Mq, 
from the A/bh — o" relation. The MGE model, constructed from 
a finite spatial resolution HST WFPC2 image, has by construction 
an unrealistic flat asymptotic density profile well inside the central 
observed pixel (r < 0'.'02) and therefore we assumed a cusp with a 
power-law slope of 1.5 (p = r~^'^) inside that radius, following 
the prescription of Emsellem et al. 1 1999). The final two-integral 
DF was computed on a fine adaptive grid (i.e. more points in the 
region of strongly-changing DF with E/Emax < 2, where Emax 
is the value of the central potential of the model excluding the black 
hole) of 140 X 79 points in {E, L^). A fine grid of LOSVDs was 
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Figure 14. Three-integral model grids of inclination angle versus M/L 
ratio, T. Contours present constant Ax^, measuring the goodness of fit 
of the models. Every symbol corresponds to a three-integral axisymmetric 
model with given inclination angle and T. The grey symbols indicate the 
best fitting models at a given inclination. The top grid presents models using 
the Full field set of kinematic constraints. The global minimum is for the 
model with i = 60° and T = 4.6. The bottom panel presents models using 
the SAURON field set of kinematic constraints. The global minimum here is 
for the model with i = 65° and T = 4.4. 



Table 6. The properties of two-integral models and comparison to three- 
integral best-fit results. 





21 model 


Full field 


SAURON field 




(1) 


(2) 


(3) 


T 


4.6 


4.6 ±0.1 


4.4 ±0.1 


i 


60° 


60° ± 5° 


65° ± 5° 



Notes - Col.(l): parameters of the two-integral model of NGC 2974; 
Col. (2): recovered parameters using the Full field spatial coverage (100" X 
100" effectively) of kinematic constraints; Col. (3): recovered parameters 
using the SAURON field spatial coverage (40" X 40" effectively) of kine- 
matic constraints. 

computed from this DF. These LOSVDs were used to compute the 
observable LOSVDs on 3721 positions, two-dimensionally cover- 
ing one quadrant of the sky plane (50" x 50"), accounting for the 
instrumental set up (size of SAURON pixels) and atmospheric see- 
ing (which matched the observations of NGC 2974 and was used 
for the three-integral models in SectionU2j. The parameters of the 
two-integral model of NGC 2974 are listed in the first column of 
Table|6l 

In order to mimic the real observations (as well as to re- 
duce the number of observables in the fit) we adaptively binned 
the spatial apertures using the Voronoi tessellation method of 



ICappellari & CopinI <2003l) as it was done for the observations, as- 
suming Poissonian noise. The final LOSVDs were used to calculate 
the kinematic moments (V, a, hs to he) by fitting a Gauss-Hermite 
series (first row on Fig. I15> . These values were adopted as kine- 
matic observables for the three-integral models. From these data 
we selected two sets of kinematic observables: 

• The first set consisted of all 513 spatial bins provided by the 
two-integral model. In terms of radial coverage this Full field set 
extended somewhat beyond two effective radii for NGC 2974. 

• The second set had a limited spatial coverage. It was limited 
by the extent of the SAURON observations of NGC 2974, approxi- 
mately covering one effective radius. We call this observational set 
of 313 bins the SAURON field. 

5.2 Recovery of input parameters 

Our two-integral models are axisymmetric by construction and all 
necessary information is given in one quadrant bounded by the 
symmetry axes (major and minor photometric axes). Hence, the in- 
puts to the three-integral code covered only one quadrant of the 
galaxy. Although only one quadrant was used for the calculations 
we show all maps unfolded for presentation purposes. 

The three-integral models were constructed in the same way 
as described in Section IO For both (Full field and SAURON field) 
sets of kinematic data we created orbit libraries of 2 x 41 x 10 x 10 
orbits and constructed models on grids of (T,i). As for the real 
data, we used a regularisation parameter A = 4. The resulting grids 
are presented in Fig. 1141 and the best-fitting models are listed in 
the Table|6| The three-integral models were able to recover the true 
input parameters within the estimated errors, although the SAURON 
field models were less accurate. 

The kinematic observables computed from the two-integral 
models are noiseless, without errors nor intrinsic scatter typical of 
real measurements. For each kinematic observable we assigned a 
constant error, but representative to the SAURON observations of 
NGC 2974: la errors for V, a, /13, and . . .he were 4 km s~^, 
7 kms~^, 0.03 and 0.04 respectively. 

As we wanted to test an ideal situation we computed the 
kinematic observables from the two-integral model without adding 
noise (intrinsic scatter) typical for real measurements. However, 
given the fact that the input model is noiseless, the levels com- 
puted from the fit to the kinematics are meaningless. In order to 
have an estimate of the uncertainties in the recovery of the parame- 
ters we computed half a dozen Monte Carlo realisations of the kine- 
matic data, introducing the intrinsic scatter to the noiseless data. For 
each Monte Carlo data sets we constructed a three-integral model 
grid like in Fig.^| Due to time limitations we calculated parameter 
grids (T, i) of models with smaller orbit libraries (2 x 21 x 7 x 7 
orbits) of both Full field and SAURON field data sets. In this case 
we also applied the regularisation scheme with A = 4. Approx- 
imate 3(T confidence levels assigned to the best-fitting parameters 
are listed in Table|6| Setting the regularisation to zero, we observed 
similar trends. 

The numbers in Table|S|suggest that the inclination is formally 
recovered by three-integral models, as seen in the case of the obser- 
vation (Section f4.3t . We repeat the exercises of plotting a sequence 
of models with different inclinations (using the Full field kinemat- 
ics). They are presented in Fig. lI6l and we can see a very similar 
trend as in Fig. 1121 there appears to be little difference between the 
models, although by scrutinising the details, it is possible to choose 
the best model by eye. 
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Figure 15. Comparison between two-integral analytical model kinematics and three-integral models. First row: two-integral model used as input to the three- 
integral code. Second row: best fitting three-integral model (i = 60° , T = 4.6) using the Full field set of kinematic constraints. Third row: best-fitting 
three-integral model (i = 65°, T = 4.4) using the SAURON field set of kinematics. The spatial extent of this set is marked by red squares on the maps. From 
left to right each panel presents: mean velocity V, velocity dispersion a, and Gauss-Hermite moments: /13, h4, hs and hg. Isophotal contours of total light are 
shown with solid lines. 



The smoothness of the data of the two-integral test model 
helps in recognising the best-fit model. The models with lower in- 
clination (towards face-on) are generally smoother than the higher 
inclination (towards edge-on) models, which systematically show 
radial structures. These "rays" visible on Fig. 1161 are artifacts of 
the discreteness of our orbit library. The starting points of the or- 
bits correspond to the positions of orbital cusps, which carry the 
biggest contribution to the observables. The total number of cusps 
is determined by the number of orbits, and the finiteness of the orbit 
library is reflected in the discrete contribu tions of the cusps to th e 
reconstructed observables (see Fig. 1.2 of ICappellari et al]|2004) . 
The projection effects, however, increase the smoothness by spa- 
tially overlapping different cusps; hence, models projected at e.g. 
i = 60° will be smoother than models viewed edge on. We believe 
this effect could influence the favouring the lower inclination 
models. 

Although present, this effect does not provide the only con- 
straint on the inclination. When examined more closely, models 
with low inclination do reproduce certain features better, for ex- 
ample: the shape and amplitude of the velocity dispersion in the 
central 20", /13 and /14 at the larger radii (towards the edge of 
the field), the central 10" of hs and h(,. In all cases the model 
with i = 60° reproduces these features better than other models. 
Comparing the models, the most significant contribution to the 
comes from the velocity dispersion, but generally, individual ob- 
servables have slightly different values, which increase with the 
inclination moving away from 60° and is visible only as a cumula- 
tive effect. This explains the similarities of the different models to 
the eye, although they are formally significantly different. 

The difference in the model observables (which include mo- 
ments up to ha) are below the level of the expected systematics 



in the data (e.g. template mismatch), or in the models (e.g. reg- 
ularisation or variations in the sampling of observables with or- 
bits). In the case of NGC 2974, using high signal-to-noise two- 
dimensional data, the difference between the models themselves 
are smaller than between the best-fitting model and the data, im- 
plying that the inclination is only weakly constrained. This result 
suggests a fundamental degeneracy for the determination of incli- 
nation with three-integral models, which is contrary to indications 
from previous work by^^rolme et al. 1 2002). Theoretical work and 
more general tests on other galaxies are needed for a better under- 
standing of this issue. 



5.3 Effect of the field coverage on orbital distribution 

The next step is to compare in more detail the three-integral models 
using the two different kinematic data sets. The kinematic struc- 
tures of the best-fitting models are presented in Fig. 1151 In the 
region constrained by the kinematic data both models reproduce 
equally well the input kinematics. As expected the regions outside 
the SAURON field are not reproduced well. It is more interesting to 
compare the phase spaces of the models. In particular, we wish to 
see whether the mass weights assigned to the orbits (represented by 
the integrals of motion) are the same for the two models. 

The corresponding integral-spaces are shown in Fig. 1191 Red 
rectangular boxes represent the extent of the kinematics used to 
constrain the models. In the regions constrained by both kinematic 
sets, the two integral spaces are identical: both models recover 
the indistinguishable orbital mass weights. The differences appear 
at larger radii (beyond 20"), outside the area constrained by the 
SAURON field kinematic set. Putting this result in the perspective 
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Figure 16. Sequence of three-integral models for different inclinations fitting the (FiiU field of kinematic observables. First row: two-integral test model. 
Subsequent rows: three-integral models for i = (55°, 60°, 70°, 80°, 89°). Columns from left to right present moments of the LOSVD: v, a, and from /13 to 
he - Isophotal contours of total light are shown with solid lines. 



of observations, the resulting phase space (in the region constrained 
by the observations) does not depend on the extent of the radial cov- 
erage used to constrain the model. This result strengthens the case 
of the NGC 2974 modelling results, where we have integral-field 
observations reaching ~ 1 re . The recovered integral space and its 
features would not change significantly if we had a spatially larger 
observational field. 



5.4 Recovery of the internal moments 

We wish to see if the best-fitting three-integral model to the 
f{E,Lz) test galaxy model is consistent with the input, i.e., 
whether three-integral models will recognise the true structure of 
the test galaxy. A first estimate can be achieved by investigating 
the internal structure of the resulting model galaxy, specifically the 
shape of the velocity ellipsoid. We define the tangential dispersion 
as at = [^i^e + ''"^)]^^^- Note that cr^ includes only random mo- 
tion so that for an isotropic distribution, under the given defini- 
tion, the radial (err) and tangential dispersion are equal. Since two- 
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Figure 17. Comparison of the integral spaces of the f{E, Lz) test models. Upper five rows belong to the model constrained by the SAURON field kinematics 
and the lower five rows to the model constrained by the Full field kinematics. The meaning of each panel is the same as in Fig. ll3l The colour bar on the right 
represents the relative colour coding of the orbital mass weights (7), where 1 represent the largest mass weight assigned to an orbit at the given energy. In the 
region constrained by both kinematic sets the integral spaces are indistinguishable. 



integral models are isotropic in the meridional plane per definition, 
we expect to recover that ar is equal to ag and the cross-term aRs 
is equal to zero. In Fig.^|we show the ratio of the moments of the 
velocity ellipsoids at different positions in the meridional plane and 
at different radii from the model constrained by the Full field kine- 
matics. One can see that = o"e within « 5%, confirming that 
our three-integral model recovers the true internal moments. Also, 
computing the cross-terms ane, we verified that it is negligible ev- 
erywhere. Similar results are also recovered from the SAURON field 



model: in the constrained region the ratio of ar and ag moments are 
consistent with unity. 



5.5 Recovery of the distribution function 

The previous result shows that the constructed three-integral model 
is a consistent representation of the input f{E, L^) model. A more 
conclusive test, however, is to compare the distribution functions. 

The results of the three-integral Schwarzschild method are or- 
bital mass weights, 73/ [E, Lz for each set of integrals of mo- 
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Figure 18. Moments of the velocity ellipsoid for the two-integral test 
galaxy recovered by the three-integral model (constrained by the Full field 
kinematic set). The upper panel shows the ratio of the radial and tangential 
moment as defined in the text. The lower panel presents the ratio of the ra- 
dial and latitudinal moments. Different lines show the ratio of the moments 
in the meridional plane at different position angles, starting from the major- 
axis to the minor axis. In a true two-integral galaxy, all three moments of 
the velocity ellipsoid have to be equal. Note that the deviation inside l"are 
expected since the data do not constrain the model in that region. 



tion {E, Lz, I3), which define an orbit. The DF is related to the 
ma ss weights via the p hase-space volume (for a detailed treatment 
see lVandervooillig'i^ 



j{E,L„l3) dEdLAh = 

f{E,L,,h) AV{E,L,,h) dEdL.dh 

where 



AV{E,L„h)^ / J{x,E,L„h)d-'x, 



(6) 



(7) 



and J{x, E, Lz, I3) is the Jacobian of the coordinate transforma- 
tion from {x, v) to {x, E, Lz, I3), and Q is the configuration space 
accessible to an orbit defined by the integrals (E, Lz, Unfor- 
tunately, J3 is not known analytically, so the above relation can not 
be explicitly evaluated, except for separable models. For this rea- 
son we limit ourselves to test the consistency of our three-integral 
mass weights with the input two-integral DF. This is possible, since 
if the recovered tiu'ee-integral DF is equal to the input DF, then the 
mass 727 (is, Lz) assigned to the stars in a given range of (E, Lz) 




Figure 19. Comparison of the mass weights from the two-integral model 
(upper panel) and the results of the three-integral modelling (lower panel) 
using the two-integral model as input. The galaxy centre is on the left side. 
The first bin containing the contribution of the black hole was not plotted, 
since the resolution of the models do not allow for its recovery. The two ver- 
tical lines enclose the region constrained by the kinematic data. The colour 
bar represents the colour coding of the orbital mass weights (7) in log scale 
assigned to any interval of (E, Lz). There are two main differences between 
the input and the best fit model: (i) the horizontal stripes, which are due to 
the discrete way in which we select orbits in Lz range, and (ii) the mass 
weights at high Lz values, which, however, contain a very small fraction of 
the galaxy mass and are not significant. 



by the three-integral model has to be equal to the mass in the same 
range of the input model. 

There exists a precise relation between the input two-integral 
DF, / — f{E, Lz), and the corresponding orbital mass weights, 
72/(£', Lz). The total mass of th e system is the integral of the DF 
over the phase-space. Using this. lCretton et aL (.1999 ) derived the 
expression for the mass weights in Appendix B of their paper (eq. 
B4): 



l2i{E,Lz 

with 

dM 
dEdL, 



dM 
dEdLz 



■dEdLz 



= f{E,Lz 



{Rdz - zdR) 



(8) 



(9) 



ZVC(E,L^) 



where the contour integral yields the area of the zero velocity curve 
(ZVC). Before applying eqs. js} and ^ on the two-integral DF we 
rebinned it to the same grid of AE and ALz as the three-integral 
mass weights. Finally, we approximate the integral js} by multiply- 
ing the mass fraction in each grid cell by AEALz. 

For our comparison we defined as the energy intervals the set 
of energies used in the construction of the three-integral models 
(total of 41). The interval in angular momentum was defined as 
a step of 0.1 of Lz/L™""= from -1 to +1 (total of 20) for a given 
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energy. The resulting grid of orbital weights is relatively coarse, but 
is representative of the model. The agreement between the two sets 
of mass weights is shown in Fig. ll9l The main features of the given 
two-integral test model are well reproduced by the three-integral 
model. Again, the mass weights should be compared in the region 
constrained by the data (between the vertical lines in the figure). In 
order to quantify the agreement we restricted the comparison to the 
part of the distribution of mass weights contributing significantly to 
the model inside this region, roughly bounded by the green levels 
on Fig. 1 191 This selected region contains about 97% of mass. The 
mean absolute deviation between the two sets of mass weights is 
6%, with peaks around 15%, except for a narrow region towards 
the edge of the kinematic coverage around E = —8 x 10'* fcm^s"^ 
and with L^/L™""' > 0, where deviations can reach 35%. The 
mass weights of the three-integral model are also relatively noisy 
which is mostly the consequence of imposed discreteness as well 
as the numerical nature of the method. 

From our test model, where the potential is known, we con- 
clude that the Schwarzschild method can reliably recover the in- 
put distribution function within the region constrained by two- 
dimensional kinematics (including higher-order moments). We 
have already shown that the existence of constraints outside of an 
effective radius does not significantly affect the orbital distribution 
inside an effective radius. This suggests that, in the region con- 
strained by integral-field kinematics, a representative DF can be 
recovered using the Schwarzschild method. 



6 MODELLING OF EMISSION-LINE GAS 

Clearly, with its prominent gas component, NGC 2974 is an un- 
usual elliptical galaxy. The observations indicate the morpholog- 
ical similarity between the small (HII) and large scale (HI) gas 
discs. Also, long-slit measurements of stellar and gas motions de- 
tect similarities betw een the stellar and gas kinematics (Kim 1989; 
ICinzano & van derMareL.1994.) . In this section we investigate the 
inclination of the gaseous component as well as construct dynam- 
ical Je ans model of the gas disc as is ICinzano & van der Marell 
^993), in order to compare them to the results of the stellar dy- 
namical modelling. 



6.1 Inclination of the gas disc 

The existence of the emission-line gas disc in NGC 2974 can be 
used to infer the inclination of the galaxy assuming an equilibrium 
dynamical configuration. This has been attempted before and in all 
studies the inclinatio n of the gas disc was con sistent with 55° — 60° 
jAmico etalJll993l 55°: iBuson et"ai]ll993l 59°; CvdM94 57.5°; 
IPlana et alJll998l 60°). The high quality of the two-dimensional 
SAURON kinematics allows us to estimate the inclination of the 
emission-line gas disc more accurately. 

We assume the motion of the emission-line gas is confined to 
a thin axisymmetric disc. We neglect the deviations from axisym- 
metry discussed in Section|3| and symmetrise the gas velocity map 
in the same way as the stellar velocity map in Section l43l This ve- 
locity map is an axisymmetric representation of the observed field, 
which can be compared to an axisymmetric model of the disc ve- 
locity maps. Constructing a model axisymmetric two-dimensional 
velocity map requires only the kinematic major axis velocity pro- 
file v'jy^j. The entire map is then given by the standard projection 




i [degrees] 

Figure 20. Ax^ as a function of inclination obtained by comparing the 
symmetrised data for NGC 2974 and the model gas disc velocity map de- 
scribed in l6. II 

formula: 

VLOs{x' ,y') ^V4.(^ ^ ^^^^ ^ ^'"mj(^), (10) 

where = x'^ + (j/'/ cosi)^ and i is the inclination of the disc. 
It is clear from this formula that for a given observed major-axis 
velocity, the velocity field is just a function of inclination. From 
the first three odd coefficients (ci, C3 and C5) of the kinemetric 
expansion we construct the velocity profile v'^j along the major 
axis {9 = (j>i). Using the major axis velocity profile, we created a set 
of disc velocity fields inclined at different values of i, and compared 
them with the symmetrised velocity field. We did not correct for 
the influence of the PSF as this effect is small and is confined to the 
central few arcseconds, which we excluded from the comparison. 
We also compared the models with the non-symmetrised velocity 
map and the results were in very good agreement, but with a slightly 
larger uncertainty range. 

Figure |20| presents the Ax^ obtained by subtracting the disc 
model velocity field from the symmetrised measurements. The 
best-fitting inclination is i = 60° ± 3° (at la level). Fig. l21l shows 
a comparison between the symmetrised and model velocity fields 
for a few representative inclinations. The differences between the 
model fields are mostly in the opening angle of the iso-velocity con- 
tours, which change with the inclination of the field. The opening 
angle of the model with i = 60° is the most similar the observed 
velocity field and this significantly lowers the of the fit. 

The best-fit inclination of 60° for the emission-line gas is in 
excellent agreement with literature values determined from the var- 
ious gas components. Our best-fitting three-integral stellar dynam- 
ical model was obtained for an inclination of 65° with Sa uncer- 
tainty of 2.5°. This inclination is close to the inclination of 60° 
presented in this section, suggesting a good agreement between the 
stellar and gaseous models. There is, however, the concern that the 
agreement may not be as significant as it seems in light of the tests 
and results from Section l5!2l 
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Figure 21. Comparisons between observed and model gas disc velocity fields. From left to right: observations (symmetrised), model fields for inclinations of 
50°, 60° (best-fitting), 70°. Isophotal contours of total light are presented with ellipses. 



6.2 A simple dynamical model for the disc 



s 



At large scale, the gas kinematic maps are consistent with the as- 
sumption that the emission-line gas is moving in a thin disc. This 
assumption clearly breaks down in the inner few arcseconds, but at 
this point we neglect this effect. The observed gas velocity disper- 
sion is high everywhere in the disc and is much larger than the ther- 
mal velocity dispersion which s hould be of the order of (Jthermai ~ 
10 kms"^ JOsterbrockl (19891) and in any case < 30 kms 
Clearly, in addition to the thermal dispersion, the gas has an- 
other source of motion, which is not presently understood, but is 
seen in many galaxies l^ Bertola et al...l995.) . Several studies (e. 
Ivan der Marel & van den Boschlll99a IVerdoes Kleiin et ai]l200l 
assumed that the non-thermal gas velocity dispersion is the re- 
sult of 'local turbulence', without describing the details of the un- 
derlying physical processes. In this assumption, the gas still ro- 
tates at the circular velocity and the invoked turbulence does not 
disturb the bulk flow of the gas on circular orbits. The alterna- 
tive to this assumption is that the non-thermal velocity disper- 
sion component comes from coUisionless gravitational motion of 
the gas, where the gas acts like stars: clumps of gas move on 
self-intersecting orbits. This is, perhaps, not very physical for gas 
in general, but it can be applied to estimate the difference be- 
tween the circular and streaming velocity (and including the pro- 
jection effect on the observed velocity) of the gas. Several stud- 
ies used this approach successfully (e.g.|c inzano & van der Marel 
1 1994 'Cretton et al."2000'; 'B arth et alj|200k .Aguerri et aL.2003^ 
iDeba ttista & Williams 2004). Presently, the role and importance of 
the asymmetric drift remains an unresolved issue. 

In constructing our simple disc model we assumed that the 
emission-line gas is moving in individual clumps that interact only 
collisionlessly. The clumps move along ballistic trajectories in a 
thin disc, under the influence of the galaxy potential given by the 
stellar distribution (Section|4j. The gas kinematics are determined 
by solvin g the Jeans equations for ra dial hydrostatic equilibrium. 
FoUowing lBinnev & Tremaind ll987l eq. 4-33), the streaming ve- 
locity of gas can be written in cylindrical coordinates as: 
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d In p d In (jj; 



R- 



dR 
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(11) 



where we have assumed the distribution function depends on the 
two classical integrals of motion, / — f{E,Lz), which implies 
an = and vrv^ = 0. In eq. Jllk Vc is the circular velocity 
(y^ -R(d$/d_R)), $ is the total potential of the galaxy obtained from 
the MGE fit assuming inclination i, gr and are the radial and 
azimuthal velocity dispersions, and p{R) is the spatial number den- 
sity of gas clouds in the disc. Lacking any alternative, we use the 




Figure 22. Fit to the surface-density profile of the [OIII] emission lines. 
Dashed lines present individual exponentials given by ea.lTll. 
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Figure 23. Ratio cru/Vc for the best-fitting model to the emission-line gas 
data. 
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surface brightness of the gas to estimate p(R). Instead of using the 
actual measured values of the [OIII] and H/3 flux, we parametrise 
the emission-line surface brightness with a double exponential law, 



p = poe "-0 + p^e 



(12) 



in order to decrease the noise. The parameters are obtained from the 
fit to the [OIII] data shown in (Fig.|22j, where the two-dimensional 
surface brightness was collapsed to a profile by averaging along 
ellipses of constant ellipticity (ellipticity of the galaxy) and position 
angle (PA of the galaxy). The errors are standard deviations of the 
measurements along each ellipse. 

The relation between radial and azimuthal velocity dispersions 
can be obtained using the epicyclic approximation. This gives (fol- 
lowing eq. 4-52 in Binne v & Tremaine 1987.) : 



1 + 



dlnK 
dlni? 



(13) 



This approximation is valid for small values of the asymmetric 
drift — Vc, or, in other words, in the limit of a cold disc with 
small velocity dispersion, cr ^ Vc. This is marginally the case in 
NGC 2974, clearly violated in the central < 5", but it is acceptable 
in most of the observed regions (Fig. l23> . 

The observed quantities can be obtained from the calculated 
intrinsic properties projecting at an inclination angle i. The pro- 
jected two-dimensional line-of-sight (LOS) velocity field is given 
by eq. <10> . Within these assumptions of the disc model, the pro- 
jected LOS velocity dispersion is: 



2 / 2 



2 ^ 



I 2 



(14) 



We constructed the asymmetric drift models of the emission- 
line gas in NGC 2974 using the MGE parametrisation of the po- 
tential, T=4.5, inclination i — 60°, and simultaneously accounting 
for the atmospheric seeing and pixel size of the SAURON observa- 
tions (see Oian et al. 1995 for details). In the process we assumed 
an exponential law for the radial velocity dispersion: 



<jR = ao + (Tie 



(15) 



Our models, therefore, have three free parameters, ao, ui, and 
Rcr, and varying them we constructed gaseous disc models of 
NGC 2974. The models were compared with the symmetrised ve- 
locity and velocity dispersion maps averaging the four symmetric 
positions on the maps as required for an axisymmetric map (Sec- 
tion UsJ. The best-fitting model was obtained for ao = 85 ± 
5 kms"\ 0-1 = 180 ± 10 kms"\ and ^ 5 ± 1". Compar- 
ison of this model with the symmetrised observations is presented 
in Fig. |24| 

This simple asymmetric drift model can reproduce rather well 
the general properties of the emission-line gas disc, including the 
bulk of the streaming velocity as well as the significant non-zero 
velocity dispersion. The overall fit is quite good with a mean dif- 
ference of only 5%, with peaks up to about 20%. The main discrep- 
ancy occurs along the major-axis at a radius around 10", close to 
where we observe the presumed elliptical ring in the ionised gas 
equivalent width (see Fig. 6). The observed minor-axis elongation 
of the gas velocity dispersion is also not reproduced by the model. 
This is not surprising since some of those features, as mentioned, 
are signatures of non-axisymmetry in NGC 2974, and cannot be 
represented by simple axisymmetric models. 

Since the normalisation of the kinematics is mostly controlled 
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Figure 24. Bottom Panels: Asymmetric drift model for the best fit- 
ting pai'ameters compared to Top Panels: the symmetrised (mirror-(anti)- 
symmetric filtering with 6 terms) mean velocity (first column) and velocity 
dispersion (second column). Overplotted isophotes are levels of the recon- 
structed total intensity. 



by the assumed mass-to-light ratio, the goodness of the fit im- 
plies that the overall energy budget coming from the Schwarzschild 
modelling is consistent with the observed gas kinematics: the high 
value of the gas velocity dispersion is therefore globally compen- 
sated by a lower mean ve l ocity. A similar finding was reported by 
ICinzano & van der MarelHl994h . All this therefore raises the ques- 
tion of how the gas can remain dynamically hot in a disc galaxy 
like NGC 2974 at a scale of ~ 1 kpc. 



7 CONCLUDING REMARKS 

This paper presents a case study of the early-type galaxy 
NGC 2974, which was observed in the course of the SAURON sur- 
vey of nearby E/SO galaxies. 

Kinematic position angles of the stellar and gaseous kinemat- 
ics of this galaxy are on average well aligned. The stellar kinematic 
maps exhibit point- and mirror-(anti)-symmetry, with the kinematic 
angle equal to the photometric PA, and are consistent with an ax- 
isymmetric intrinsic shape. The gaseous velocity map is more com- 
plicated, with clear departures from axisymmetry in the centre of 
the galaxy (< 4"). At larger radii, the gas kinematic angle is not 
constant, although is largely consistent with the photometric PA, 
showing deviations of a few degrees. The SAURON observations 
of NGC 2974 confirm the existence of non-axisymmetric pertur- 
bations consistent with an inner bar (EGF03) as well as a possible 
large-scale bar. The departures from axisymmetry are not visible in 
the stellar kinematic maps and therefore are likely to be weak. This 
allows us to construct axisymmetric models of NGC 2974. 

We constructed self-consistent three-integral axisymmetric 
models based on the Schwarzschild' s orbit superposition method, 
varying mass-to-light ratio and inclination. The observed surface 
brightness was parameterised by multi-Gaussian expansion model 
on both ground- and space-based imaging. The models were com- 
pared with the SAURON kinematic maps of the first six moments 
of the LOSVD (v, a, hs - ha). The best-fitting model has T — 
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4.5 ±0.1 and i = 65° ± 2.5°. The inclination is formally well 
constrained, but there are several indications that the recovery of 
the inclination is uncertain: (i) differences between the models are 
on the level of the expected systematics in the data (e.g. template 
mismatch); (ii) difference between the best-fit model and the data 
are bigger than differences between other models and the best-fit 
model; (iii) limitations of the models (discreteness of the cusps) 
may artificially constrain the inclination. 

The internal structure of NGC 2974 (assuming axisymmetry) 
reveals the existence of a rapidly rotating component contributing 
with about 10% of the total light. This component is composed 
of orbits allowing the third integral and does not represent a cold 
stellar disc, although suggests a flattened structure similar to an SO 
galaxy. 

The results of the stellar dynamical models were compared 
with the results of modelling the gas component. The inclination 
of the gas disc, calculated from the emission-line velocity map is 
i = 60° ± 3, in agreement with the formally constrained stellar 
inclination. A simple model of the gas disc in the same potential 
used for the stellar modelling (T = 4.5, but i — 60°) was able 
to reproduce the characteristics of the gas kinematics (v and cr) on 
the large scale, but, as expected, failed in the centre. In general, 
the observed gas kinematics is consistent with being produced by 
gravitational potential determined by the stellar dynamical models, 
but the origin of the high velocity dispersion of the gas on large 
scales remains an open question. 

We performed a set of tests of our implementation of the 
Schwarzschild's orbit superposition method. For this purpose we 
constructed a general two-integral model of NGC 2974, and used 
the reconstructed kinematics as inputs to the Schwarzschild's 
method. We tested (i) the influence of the radial coverage of the 
kinematic data on the internal structure, (ii) the recovery of the test 
model parameters (T,i), and (iii) the recovery of the test model DF. 
The tests show that: 

(i) Increasing the radial coverage of the kinematic data from Ire 
to 2re does not change the internal structure within Ire. The results 
of the dynamical models of the S AURON observations of NGC 2974 
would not change if the radial coverage would be increased by a 
factor of 2. 

(ii) We find that three-integral models can accurately recover the 
mass-to-light ratio. Although the models are also able to constrain 
the inclination of the test model formally, the apparent differences 
between the models are small (as in the case of the real observa- 
tions). Under careful examination, it is possible to choose the best 
model by eye, but the decisive kinematic features are below (or 
at the level) of the expected systematics in the data (e.g. template 
mismatch) and might be influenced by the uncertainties in the mod- 
els (e.g. regularisation or variations in the sampling of observables 
with orbits). This suggest a degeneracy of models with respect to 
the recovery of inclination. More general tests on other galaxies and 
theoretical work is needed for a better understanding of this issue. 

(iii) From a realistic test model, the analytically known input DF 
is recovered within 6% in the region constrained by integral-field 
kinematics. This suggests that applying the Schwarzschild tech- 
nique with integral-field kinematics can reliably recover a repre- 
sentative DF. 
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